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Abstract: A study was conducted to examine the distribution pattern of 
four rhizomatous medicinal and aromatic plant species (MAPs) viz., 
Angelica glauca, Pleurospermum angelicoides, Rheum emodi and Arne- 
bia benthamii in different forest stands in Central Himalaya. Results 
show that A. glauca and P. angelicoides had a higher (50%) frequency at 
Chipkoan, Garpak and Phagati forest, R. emodi had a higher (60%) fre¬ 
quency at Rishikund, Suki and Himtoli, and A. benthamii had a higher 
(70%) frequency at Suki and Khambdhar The densities of A. glauca (0.6 
plants-m' 2 ) and P. angelicoides (0.5 plants-m' 2 ) were higher at Chipkoan 
and Garpak sites than at other micro-sites, while densities of R. emodi 
(0.8 plants-m' 2 ) and A. benthamii (1.0 plants-m' 2 ) were higher at Suki and 
Khambdhar sites. A. glauca had highest total basal covers (TBC) (1.2 
cm 2 -m' 2 ) at Chipkoan, P. angelicoides had highest TBC (0.92 cm 2 -m' 2) at 
Lati kharak site, A. benthamii had the highest TBC (6.48 cm -m') at 
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Khambdhar, and R. emodi had highest TBC (4.53 cm 2 -m' 2 ) at Rishikund. 
For the four studied species, A. glauca showed a contagious distribution, 
P. angelicoides and R. emodi showed the random and A. benthamii 
showed the regular type of distribution. 

Keywords: alpine ecosystem; Himalaya; medicinal and aromatic plants; 
traditional knowledge; Uttarakhand 

Introduction 

Among many biological hotspots around the world, the Hima¬ 
laya and Western Ghats in the Indian sub-continent are the bio¬ 
diversity rich areas. Himalaya covers only 18% of the Indian 
sub-continents; however, it accounts for more than 50% of In¬ 
dian forest, of which 40% are endemic species (Maikhuri et al. 
2000; Kala 2005). It is well known that the distribution pattern of 
plant species is mainly governed and regulated by altitude as 
well as edaphic and climatic factors (Bongers et al. 1999; Nauti- 
yal et al. 2001; Kala 2004; Kharkwal et al. 2005, 2007) and their 
population is especially affected by human activities due to 
various reasons (Niggemann et al. 2009). Himalayan forests 
extend from the lower mountain slopes with nearly tropical to 
temperate timberline (Singh and Singh 1992). The alpine region 
of the Himalaya has a variety of medicinal plant species. How¬ 
ever, exploitation through legal or illegal means has ultimately 
resulted in the decline of the population of many valuable me¬ 
dicinal and aromatic plants (MAPs) that are of ecological and 
economic significance. 

Moreover, a large number of medicinal plants have become 
threatened due to their small population size, narrow distribution 
area, habitat specificity, destructive mode of harvesting, heavy 
livestock grazing, high value of utilization, climate change, ha¬ 
bitat loss, present development activities and genetic drift (Kala 
2000; 2005). Although extensive research has been carried out 
on many medicinal plants of this region (Samant 1993; Purohit et 
al. 2001; Kumar and Ram 2005; Misra et al. 2008; Kala 2006, 
2010), there is relatively less information available about the 
distribution pattern of large number of MAPs. In this regard, the 
present study was initiated to investigate the distribution pattern 
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and conservation measures of four rhizomatous species {Angel¬ 
ica glauca Edgew., Pleurospermum angelicoides (DC.) C.B. 
Clarke, Rheum emodi Wall, ex Meissn., and Arnebia benthamii 
Wall, ex G. Don.) in the Dhauli Ganga catchment of Nanda Devi 
Biosphere Reserve in Central Himalayan region of India. 


and the soils are loam to sandy loam. The study area is domi¬ 
nated by the Bhotiya; a tribal community belonging to an In- 
do-Mongoloid ethnic group (Fig. 1). 

Climate 


Materials and methods 

Study area 

The study was conducted in Dhauli Ganga catchment (30° 16' to 
30°32' N and 79°44' to 80°02' E) of Nanda Devi Biosphere Re¬ 
serve (NDBR), which lies in Chamoli district of Uttarakhand 
(Central Himalaya). The area is dominated by crystalline rocks 


The daily minimum and maximum temperatures were recorded 
from March 2000 to March 2003. Monthly minimum tempera¬ 
ture ranged from 2.2°C to 16°C, and the monthly maxi¬ 
mum-temperature ranged from 15.3°C to 27.2°C. June and Au¬ 
gust were the warmest months of the year, with an average tem¬ 
perature of 27°C. Average rainfall is about 937 mm/year and 
about 43% of annual rainfall occurs over a short period of two 
months (July-August), with a strong monsoon influence (Fig. 2). 



Fig. 1 Location map of the study site in the Nanda Devi Biosphere reserve (Not to Scale), India 


Survey 

Specimens were collected with the help of ‘Flora of Chamoli’ 
(Naithani 1984) and information on the botanical work of the 
Nanda Devi National Park (Hajra, 1983). Plant species were 
identified at the herbarium of Botany Department of HNB 
Garhwal University. Plant specimens were collected, preserved 
and identified following the techniques reported by Gaur (1999) 
and Naithani (1984) and specimens were deposited in the her- 
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barium. 

Considering the importance of these endangered taxa in the 
Himalayan region, complete and in-depth quantitative field sur¬ 
veys were conducted on population estimation in different locali¬ 
ties. The survey was undertaken in a wide transverse range cov¬ 
ering more than 45 km in the study area (from Raini to Malari in 
Chamoli District) in Uttarakhand. Five prominent sites between 
elevations of 2 000 and 3 500 m a.s.l. were analyzed for each 
species (Table 1). Approximately 100 m x 100 m forest stands 
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were selected in the area. In each forest stand, 10 quadrats (1 m x 
1 m size each) were laid randomly following the methodology of 
Kershaw (1973). Analytical features of the population such as 
percent frequency, density and total basal cover were calculated 
using the standard methods of Misra (1968). The ratio of abun¬ 
dance to frequency (A/F), a relative measure to present the dis¬ 
tribution of species in a community, was calculated according to 
the methods of Curtis and Cottom (1956) as: A/F <0.025 (regu¬ 
lar), between 0.025 and 0.05 (random), and > 0.05 (contagious) 
distribution. Importance Value Index (IVI) of each species was 
calculated for the determination of dominance and ecological 
success of a species (Curtis and McIntosh 1950). The plant spe¬ 
cies diversity was determined using Shannon-Weiner informa¬ 
tion index (Shannon and Weaver 1963). The community similar¬ 
ity coefficient of the stands was calculated on the basis of IVI 
following the methods of Jaccard (1912). Niche width for each 
species was computed using the equation given by Levins (1968). 
The degree of presence (P) of a species was determined follow¬ 
ing Braun-Blanquet (1932) method. 




.03 

M-H 

Ch 

Dh 




Rainfall 

Minimum temperature 



Maximum temperature 


Fig. 2 Climatic data for the study area, based on average rainfall and 
temperature (minimum and maximum) 


Table 1. Site characteristics of the selected population 


Sites 

Altitude 

(m a.s.l.) 

Aspects 

Slope 

(°) 

Habitat 

Khambdhar 

3900 

North West 

40-50 

Rocky terrain, moist gras¬ 
sy slopes. 

Himtoli 

3700 

South West 

50-55 

Steep slopes very moist 
grassy slopes. 

Rishikund 

3900 

North West 

45-60 

Rocky slope on Northwest 
facing slopes. 

Tolma forest 

3200 

North West 

40-45 

Mixed forest (Cedrus, 
Betula spp). 

Suki 

3300 

North West 

50-55 

Dry rocky crevices 

Chipkoan 

3500 

South East 

50-60 

Moist rocky, grass domi¬ 
nated 

Garpak 

3200 

South East 

45-55 

Moist rocky grass domi¬ 
nated 

Phagati forest 

3500 

South East 

50-55 

Moist laden moist rocky, 

boulders 

Lati kharak 

3700 

South East 

50-60 

Gentle slope in moist with 

Cedrus forest 


Results 

Occurrence and availability 

The distribution of species in the study areas was found mostly 
habitat specific (Table 2). A. glauca grows well in moist and 
shady areas especially in alpine slopes. The frequency of A. 
glauca ranged from 30% at Phagati forest site to 50% in all the 
study sites. The density of A. glauca was higher (0.6 plants-m' 2 ) 
at Chipkoan and Garpak sites. The species accompanied with A. 
glauca were Veronica cana , Fragaria nubicola and 
Calamagrostis pilosula. F. nubicola was a dominant species in 
all the sites, although, other species were also present as 
co-dominant species at different sites. P. angelicoides was 
mostly found growing in shady places of the forest. Frequency of 
P. angelicoides ranged from 30% at Lati kharak site to 50% at 
Chipkoan, Garpak and Phagati forest sites. Density of P. 
angelicoides was higher (0.5 plants-m" 2 ) at Chipkoan and Garpak 
sites. Common species growing along with P. angelicoides were 
Sedum imbricatum, Thalictrum alpinum and Polygonum spp 
(Table 3). R. emodi grows well in moist rocky places between 
boulders and stream (with 40°-50° steepness); the frequency of 
occurrence of this species was 40% at Tolma forest site and 60% 
at Sulci, Himtoli and Rishikund sites. The density of R. emodi 
was the highest (0.8 plants-m' 2 ) at Suki site. The common species 
associated with R. emodi were Danthonia cachemyriana , 
Ligularia arnicoides and Potentilla fulgens. D. cachemyriana 
was found to be the most dominant species across all sites, 
although different species were present as co-dominant species at 
various sites (Table 3). For A. benthamii, frequency of 
occurrence was 20% at Tolma forest, 70% at Suki, and 
Khambdhar, and 30% at Himtoli and Rishikund sites. The 
density of A. benthamii was highest (1 plants-m" 2 ) at Khambdhar 
site. In general, species such as D. cachemyriana , P. 
polystachyum and Fritillaria roylei were dominant at all the 
investigated sites. A/F ratio presents the idea of distribution 
patterns of species in a community. A. glauca shows contagious 
distribution, P. angelicoides and R. emodi shows random and A. 
benthamii shows regular type of distribution. A. glauca exhibited 
the highest total basal cover (TBC) (1.2 cm 2 -m" 2 ) at Chipkoan, P. 
angelicoides (0.92 cm 2 -m' 2 ) at Lati kharak, R. emodi (4.53 
cm 2 -m" 2 ) at Rishikund and A. benthamii (6.48 cm 2 -m' 2 ) at 
Khambdhar sites (Table 3). 

The maximum IVI of R. emodi (25.63) , A. benthamii (65), A. 
glauca (28.02), and P. angelicoides (33.06) were found at Him¬ 
toli, Khambdhar, Garpak, and Chipkoan sites, respectively (Ta¬ 
ble 3). In the present study, the value of diversity index (H') 
ranged from 1.47 to 2.52 for all the four species (Fig. 3). The 
highest community similarity for A. glauca (81.8) and P. angeli¬ 
coides (94.7) was both at Garpak site. Similarly, the highest si¬ 
milarity coefficient of R. emodi (90.9) and A. benthamii (90) was 
observed to occur at Khambdhar site (Table 4). Percent contribu¬ 
tion to community basal cover for R. emodi ranged between 
4.568 at Suki, to 16.566 at Khambdhar, similarly A. benthamii 
had community basal cover between 3.61 at Tolma forest and 
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81.18 at Suki,, A, glauca between 0.86 at Phagati forest and 2.11 forest and 20.35 at Rishikund (Table 5). 
at Lati kharak, and P. angelicoides between 0.151 at Phagati 


Table 2. Characteristic of Angelica glauca, Pleurospermum angelicoides, Arnebia benthamii and Rheum emodi 


Species 

Habitat and occurrence 

Altitude 

(m a.s.l) 

Life span 

Flower 

Leaves 

Root 

Angelica glauca 

Moist shady, sub-alpine slopes near, water 
springs and occur mostly in the shady 
canopy of Betula utilis forest 

1800-3700 

Perennial 

Compound umbel yellowish or purple in 
colour 

1-3 pinnate 

Tuberous 

Pleurospermum. 

angelicoides 

Mostly forest; cover and shady slope. 

2800-3500 

Perennial 

Compound umbel yellowish or purple in 

colour 

Broad leaves bilat¬ 
erally divided 

Tuberous 

Rheum emodi 

Moist rocky places, between boulders and 3500-5000 

stream. 

Perennial 

Axillary panicles, dark purple or pale orbicular or broadly 

red and tall ovate 

Stout 

Arnebia benthamii 

Sub-alpine forest and shady slopes 

3000-3900 

Perennial 

Flowers are pink to purple or maroon, in 
dense, terminal spikes; nutlets ovoid. 

Narrow-lanceolate 

Stout 


Table 3. Phytosociological attributes of Angelica glauca, Pleurospermum angelicoides, Rheum emodi , and Arnebia benthamii in the study area. 


Species 



Phytosociological attributes 



Sites 

Frequency 

Density 

Abundance 

Total basal 

Importance 

Dominant associates 


F (%) 

(Plants nf 2 ) 

frequency ratio 

2 2 

cover ( cm m" ) Value Index 


Chipkoan 

50 

0.6 

0.02 

1.2 

8.61 

C. pilosula, V. cana 

Garpak 

50 

0.6 

0.60 

0.97 

28.02 

F. nubicola, Cicerbita cyanea 

Angelica glauca Phagati forest 

30 

0.3 

0.03 

0.60 

8.01 

F. nubicola, P. arbuscula 

Rishikund 

50 

0.5 

0.02 

0.64 

13 

C. pilosula, V.cana 

Lati kharak 

50 

0.2 

0.05 

0.57 

9.36 

V. cana, C. pilosula 

Chipkoan 

50 

0.5 

0.02 

0.36 

33.6 

P. spp, S. imbricatum 

Garpak 

50 

0.5 

0.02 

0.36 

9.86 

P. spp, S. imbricatum 

Pleurospermum angelicoides Phagati forest 

50 

0.3 

0.03 

0.21 

5.74 

S. imbricatum, C. pilosula 

Rishikund 

40 

0.4 

0.02 

0.02 

7.88 

S. imbricatum, L. arnicoides 

Lati kharak 

30 

0.3 

0.03 

0.92 

12.27 

S. imbricatum, T. alpinum 

Suki 

60 

0.8 

0.02 

2.88 

14.56 

D. cachemyriana, P.fulgens 

Khambdhar 

50 

0.6 

0.02 

4.33 

24.34 

D. cachemyriana, F. roylei 

Rheum emodi Himtoli 

60 

0.6 

0.01 

4.33 

25.63 

D. cachemyriana, L. arnicide 

Rishikund 

60 

0.7 

0.019 

4.53 

16.09 

D. cachemyriana, P.fulgens 

Tolma forest 

40 

0.4 

0.025 

2.88 

17.55 

D. cachemyriana, P.fulgens 

Suki 

70 

0.9 

0.01 

5.83 

18.44 

P. polystachyum, A. margaritacea 

Khambdhar 

70 

1.0 

0.02 

6.48 

65.07 

D. cachemyriana, P. polystachyum 

Arnebia benthamii Himtoli 

30 

0.4 

0.44 

2.88 

38.56 

P. astrosanguinea, F. roylei 

Rishikund 

30 

0.3 

0.03 

2.16 

14.79 

D. cachemyriana, F. roylei 

Tolma forest 

20 

0.2 

1.44 

1.44 

8.63 

L. arnicoides, D. cachemyriana 
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Fig. 3 Species diversity (H) for 1) Rheum emodi, 2) Arnebia ben¬ 
thamii, 3) Angelica glauca and 4) Pleurospermum angelicoides. 

SK—Suki, KD—Khambdhar, HT—Himtoli, TF—Tolma forest, 
CK-Chipkoan, GP—Garpak, PF—Phagati forest, RK—Rishikund, LK—Lati 
kharak 
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Niche width value based on density value ranged from 1 to 8.2 
among the different species. Niche width value was 1.68 for A. 
glauca , 3.7 for A. benthamii , 4.7 for R. emodi , and 2.2 for P. 
angelicoides , respectively. However, Percent presence (% P) of 
all the four species was very high which might be due to their 
specific natural requirements (Fig. 4). 


Discussion 

A. glauca showed contagious distribution and P. angelicoides 
and R. emodi presented the random and A. benthamii was the 
regular type of distribution, which is a common feature in alpine 
vegetation. These patterns are consistent to the previous observa¬ 
tions (Singh and Yadav 1974; Saxena and Singh 1980). Accord¬ 
ing to Odum (1971), contagious distribution is the most perva¬ 
sive pattern in nature; random distribution is confined only in 
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very uniform environments whereas regular distribution occurs 
in those areas where competition among several individuals ex¬ 
ists. Contagious distribution depends on the local habitat, daily 
and seasonal weather change and reproductive process. 



Fig. 4 Percent presence and Niche width of the species in the study 
area. 


Table 4. Community similarity coefficients between different popu¬ 
lations of Rheum emodi , Arnebia benthamii , Angelica glauca, and 
Pleurospermum angelicoides in the study area 


Stands 

Suki 

Khambdhar 

Himtoli 

Rishikund 

Tolma forest 

R. emodi 






Suki 

100 

90.9 

73.68 

60 

80 

Khambdhar 

— 

100 

73.68 

60 

80 

Himtoli 

— 

— 

100 

58.82 

70.5 

Rishikund 

— 

— 

— 

100 

55.5 

Tolma forest 

— 

— 

— 

— 

100 

Stands 

Suki 

Khambdhar 

Himtoli 

Rishikund 

Tolma forest 

A. benthamii 






Suki 

100 

90 

73.66 

73.68 

66.66 

Khambdhar 

— 

100 

82.35 

70.58 

73.68 

Himtoli 

— 

— 

100 

75 

77.77 

Rishikund 

— 

— 

— 

100 

66.66 

Tolma forest 

— 

— 

— 

— 

100 

Stands 

Chipkoan 

Garpak 

Phagati 

Rishikund 

Lati kharak 




forest 



A. glauca 






Chipkoan 

100 

81.81 

54.54 

45.45 

47.61 

Garpak 

— 

100 

55.55 

44.44 

47 

Phagati forest 

— 

— 

100 

77.77 

70.58 

Rishikund 

— 

— 

— 

100 

58.82 

Lati kharak 

— 

-v 


— 

100 

Stands 

Chipkoan 

Garpak 

Phagati 

Rishikund 

Lati kharak 




forest 



P. angelicoides 






Chipkoan 

100 

94.73 

80.00 

60.00 

82.35 

Garpak 

— 

100 

73.68 

52.63 

75.00 

Phagati forest 

— 

— 

100 

60 

70.58 

Rishikund 

— 

— 

— 

100 

70.58 

Lati kharak 

— 

— 

— 

— 

100 


Table 5. Contribution of the species community dominance/TBC a) 
Rheum emodi b) Arnebia benthamii, c) Angelica glauca and d) Pleu¬ 
rospermum angelicoides. 


Sites 

Community 

% Contribution of 

Community 

Contribution of the 


Density 

the Species to 

total 

species community 


(Plants m‘ 2 ) 

community Density 

basal cover 

dominance/TBC 

Rheum emodi 




SK 

594.8 

0.1344 

62.78 

4.5868 

KD 

272.3 

0.2203 

26.167 

16.5662 

HT 

442.0 

0.13574 

34.401 

12.586 

RK 

401 

0.1746 

87.14 

5.2048 

TF 

301.3 

0.132 

4.52 

8.5624 

Arnebia benthamii 




SK 

827 

0.1088 

71.85 

81.1899 

KD 

394.8 

0.2532 

12.232 

52.972 

HT 

310.8 

0.1287 

9.592 

30.1059 

RK 

421.7 

0.0711 

27.407 

7.9030 

TF 

309.7 

0.0645 

39.934 

3.6159 

Angelica glauca 




CK 

249.3 

0.2406 

78.6925 

1.52492 

GP 

88.7 

0.6764 

6.184 

1.5716 

PF 

166.8 

0.1798 

69.734 

0.8604 

RK 

136.9 

0.1158 

31.0382 

2.0619 

LK 

136.9 

0.146 

26.447 

2.1179 

Pleurospermum angelicoides 



CK 

158 

0.3164 

145.199 

0.2479 

GP 

146.2 

0.3419 

133.726 

0.2692 

PF 

207.8 

0.1443 

142.613 

0.151 

RK 

54.7 

0.2585 

111.919 

20.3502 

LK 

72.4 

0.4143 

80.600 

1.1445 

Note: 

SK—Suki,, 

KD—Khambdhar, 

HT—Himtoli, 

TF—Tolma forest, 


CK-Chipkoan, GP—Garpak, PF—Phagati forest, RK—Rishikund, LK—Lati 
kharak 

The higher IVI value indicates that most of the available re¬ 
sources are being utilized by that species and the residual re¬ 
sources are being trapped by another species as the competitors 
and associates (Kukshal et al. 2009). In the present study, the 
value of diversity index ranged from 1.4 to 2.5, which is within 
the range of general temperate forest from 1.16 for young stands 
to 3.4 for old stands (Knight 1975). Niche width is a measure of 
different resources used by an organism which further provides 
the role of species on a community and the degree of specializa¬ 
tion of a species as its ability to exploit an environment range in 
space (Bisht and Kusumlata 1993). It was found that D. cache- 
myriana exhibits a broader niche width as compared to other 
species. According to Smith (1980, 1990) the species with wider 
niche are considered to be more generalized. MacArthur (1965) 
pointed out the utility of niche is to determine differences be¬ 
tween individuals and species. 

Song et al. (1997) stated that the distribution and diversity of 
plant species in a landscape depend on various factors (e.g., dis¬ 
persal, ability, competition, environmental factors such as solar 
radiation, temperature and soil geological conditions) which may 
influence the landscape vegetational structure and would show 
significant effects on richness and diversity (Heydari and Mad- 
havi 2009). However, some species are widely distributed across 
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a range of elevations, 3000 m. a.s.l. Kala (2005) reported eleva¬ 
tion as a strong feature for establishment of plant communities. 
The present findings are relevant to sustainable harvesting of 
species in their natural habitats. The observed parameters on 
population status and habitat preference would assist in under¬ 
standing the ecology and development of a conservation plan 
with regard to A. glauca , A. benthamii , R. emodi and P. angeli- 
coides. Considering the above facts, besides protecting these 
plants in their natural habitats, sustainable harvesting of these 
selected MAPs is required as per the demands of the local people 
so as to improve their socio-economic conditions. 
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